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C ~ 3 The role of neutrinos in big bang nucleosynthesis is discussed. The bounds on the 

, number of neutrino families, neutrino degeneracy, parameters of neutrino oscilla- 

tions are presented. A model of chemically inhomogeneous, while energetically 
■ smooth, universe created by inhomogeneous cosmological neutrino asymmetry is 

described. Nucleosynthesis limits on production of right-handed neutrinos are con- 
sidered. 
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1 Introduction 



Big Bang Nucleosynthesis (BBN) is one of the strongest evidences in favor of the 
Standard Hot Cosmological Model. The theory predicts abundances of light ele- 
\Q , ments 2 H , 3 He, 4 He, and 7 Li which span 9 orders of magnitude in a good agreement 

with observations. Theoretical calculations arc well defined and with the existing 
uncertainties in the cross-sections of the relevant nuclear reactions the accuracy of 
calculations is about 10% or, depending upon the clement, is even better, espe- 
cially for 4 ff e. However, comparison of theoretical results with observations is not 
straightforward because the data are subject to poorly known evolutionary effects 
and systematic errors. Still, even with these uncertainties, BBN permits to elim- 
inate many modifications of the standard model and to derive strong restrictions 
on properties of elementary particles, in particular, on neutrino. 

In what follows I will briefly discuss physics of BBN and essential parameters 
that determine production of light elements (sec. 2). After that the role of neutrinos 
in BBN is described. The bound on the number of neutrino species is presented in 
section 3. In section 4 neutrino degeneracy and its impact on BBN are discussed, 
jj] ■ Neutrino oscillations and their possible influence on BBN are considered in sec. 

5. Inhomogeneous BBN related to possible spatial variations of lepton asymmetry 
and producing large fluctuations of primordial abundances in different regions of 
the universe is discussed in sec. 6. In sec. 7 limits on some other neutrino properties 
(mass, magnetic moments, right-handed currents, etc) are presented. 

The frameworks of the talk do not permit to give the complete list of references 
to the discussed problems, so I indicated only original and most recent papers. A 
history of the problems and much more detailed list of references can be found in 
the review u. 

2 Physics of BBN 

Building blocks for production of light nuclei were prepared in the universe when 
the temperature was about 1 MeV. The reactions 

n + e + «-» p + v e , 



ibbn: submitted to World Scientific on February 1, 2008 



n + v e «-» p + e (1) 

maintained thermally equilibrium ratio of neutrons to protons 

n/p = exp(-Am/T - &) (2) 

till the temperature T dropped down below Tf = 0.6 — 0.7 MeV. Here Am = 1.3 
MeV is the neutron- proton mass difference and £ e = /i e /T is the dimensionless 
chemical potential of electronic neutrinos. At temperature below Tf the reactions 
(||) became practically frozen, their rate T ~ G F T 5 became much smaller than the 
universe expansion rate, H = a/a ~ T 2 , where Gf is the Fermi coupling constant, 
H is the Hubble parameter and a(t) is the cosmological scale factor. Correspond- 
ingly the ratio n/p would remain constant if neutrons were stable. Due to neutron 
decay with the life-time r n = 888 ± 2 sec the n/p-ratio drops as exp(— t/r n ). 

Formation of light elements through the chain of reactions p + n — > 2 H + 7, 
p + 2 H -> 3 iJe, n + 2 H ^ 3 H, 3 He + n^ A He, etc started at T — T NS = 60-70 
keV. The concrete value of Tjvs depends upon baryon-to-photon ratio 77 = ub /n~/ ~ 
5 • 10 -10 (found from BBN itself and now from the angular spectrum of CMBR as 
well) and the nuclear binding energies. The temperature Tjvs is so much smaller 
than the nuclear binding energy because of very large number of cosmic photons per 
one baryon. At T — T^s all neutrons were quickly captured and no free neutrons 
remained in the plasma. Since 4 He has the largest binding energy practically all 
the neutrons, that survived to the moment when T = Tjvs, ended their lives in A He 
nuclei. 

The universe cooling rate can be obtained if one compares the cosmological 
energy density expressed in terms of the Hubble parameter and the energy density 
of thermally equilibrium plasma: 

_ 3H 2 m pl _n 2 

Ptot ~ ~^r~ - m 9 * T (3) 

where 5, counts the number of relativistic species in the primeval plasma; it includes 
2 from photons, 4 • (7/8) from e ± -pairs, and 2 ■ 3 • (7/8) from three light neutrino 
families: 

.9* = 10.75+ ? - {N v -3), (4) 

and the last term describes contribution of some non-standard energy. The latter 
could be additional neutrino species, possibly sterile, non-coupled to W and Z 
bosons, or new abundant massive particles, or some unknown form of dark energy 
(e.g. vacuum energy, or in other words, cosmological constant). In the last two 
cases equation of state of these unknown forms of matter would be different from 
relativistic equation of state p = p/3 valid for light neutrinos and thus the parameter 
N v could be a function of time. Moreover, the effect of this parameter on production 
of different elements may differ from the effect induced by additional neutrinos. 
As follows from eq. (j^) the universe cooling rate is determined by the expression: 

secy Vmcv; \ 9* ) 
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Since the cooling rate depends upon the effective number of particle species, BBN 
is sensitive to any form of energy present in the cosmic plasma in the range of 
temperatures from a few MeV down to 60 keV. 

As one can see from the discussion above primordial abundances of light ele- 
ments are the functions of: 

1. weak interaction rate, which determines the moment when the reactions of 
neutron-proton transformation freezes; this rates is determined by the neutron 
life-time; 

2. cosmological energy density parametrized as {N v — 3); 

3. number density of baryons, 7710 = 10 10 rts/n 7 ; 

4. neutrino degeneracy, given by dimensionless chemical potentials, £ a = (J. a /T, 
where a = e, fi, r; £ a remains constant during adiabatic expansion if corre- 
sponding leptonic charge is conserved. 



3 Number of neutrino families 



As we have already mentioned in the previous section, an addition of an ex- 
tra neutrino family would change the freezing temperature of rt/p-transformation 

(El). Comparing reaction rate and the rate of cosmological expansion one can find 
' — ■ 1/6 

Tf ~ 3* .So with rising number of species the number density of neutrons and 
correspondingly the mass fraction of produced A He also rises. Another effect lead- 
ing in the same direction is that the time when the nucleosynthesis temperature, 
Tjvs is reached also depends upon g„; according to eq. @ ijvs ~ .9* so with 
larger less neutrons would decay before the onset of nucleosynthesis and more 
A He would be produced. Both these effects give rise to the increase of mass fraction 
of primordial hclium-4 by approximately 5% if one extra neutrino family is added. 

The dependence of helium production on the number of neutrino families was 
first mentioned in ref. □ and a little later in ref. u. A detailed investigation was 
done in the papers tm 

Comparing observational data with the theory one can deduce an upper limit 
on the number of additional neutrino families, AJV„ = N v — 3. According to the 
numerous literature on the subject, this limit oscillates between 2 and 0.1. The 
recent analysis of ref. El gives N v = 1.8 — 3.9. 

It is usually assumed that neutrinos have the equilibrium spectrum: 



/ I/ = [exp(£/T-0 + l] 1 
and the corresponding energy density of neutrinos plus antincutrinos is 



(6) 



Pv = —-9 / dpp 
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The limit on the number of neutrino families presented above is obtained under 
assumption that neutrinos are not degenerate, i.e. £ = 0. This limit is also valid for 
non-equilibrium and v T with the same energy density as the equilibrium ones. 
This is not so for electronic neutrinos because the latter effect the n/p-ratio not only 
by their energy density but also more directly by their spectrum, since they took 
part in the reactions (]]]). Surprisingly the spectrum of all neutrinos noticeably 
deviates form the equilibrium one, though they may be exactly massless. This 
deviation is induced by different temperatures, of electrons and neutrinos due to 
e + e~-aimihilation after neutrino decoupling 013. According to analytical estimate 
of ref. Q the spectral distortion has the form: 

_.o m-4 E (HE 



- 3 • f (T 4 3 (8) 

where Sf — f — /( e? ). The distortion of the spectra of !/„ and v T is approximately 
twice weaker. The shift of helium-4 mass fraction due to this neutrino heating by 
the residual e + e~-annihilation is only a fewx 1CP 4 , though the total neutrino energy 
density becomes larger than the standard one by about 3% □. There is another effect 
of the similar magnitude and sign E3, namely finite-temperature electromagnetic 
corrections to the energy density of 7e + e~-plasma. It adds 0.01 effective number 
of extra neutrino species. 

4 Lepton Asymmetry 

If the number density of particles is different from the number density of antiparti- 
cles, i.e. charge asymmetry is non-vanishing, it is described by a non-zero chemical 
potential, £ (ph. Of course this description is valid only in the case of kinetic equi- 
librium when the distribution in energy has the canonical form dictated by kinetic 
equation with strong elastic scattering term. Fast annihilation processes imply also 
opposite values of chemical potentials for particles and antiparticles: £ = — £. 

According to eq. (||) energy density of degenerate neutrinos in thermal equi- 
librium is larger than the energy density of non-degenerate ones. The role of de- 
generate and v r in BBN is simply to increase the total energy density which 
corresponds to 



a—fi.T 



£ \ 4 ft x 2 



(9) 



The bound on AA„ presented above can be translated to the bound on the magni- 
tude of the chemical potentials. In particular, if AA^ < 1 then |£^, T | < 1.5. 

The bound on the value of |£ e | is much stronger because degeneracy of elec- 
tronic neutrinos exponentially shifts n/p-ratio, see eq. (§). The BBN bounds on 
chemical potentials would be somewhat weaker if combined variation of all chem- 
ical potentials is allowed. In this case_a large value of \^, T \ may be compensated 
by a relatively small and positive £ e til. Recent analysis of the work 113 based on 
additional information extracted from the measurements of angular fluctuations of 
CMBR permits to obtain the limits: 

-0.01<£„ o <0.2, |6, M ,„ T |<2.6 (10) 
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under assumptions that the primordial fraction of deuterium is_.D/_ff = (3.0 ±0.4) • 
10~ 5 . More details and references can be found in the review 111. 

These results can be further strengthen if there are neutrino oscillations that mix 
24, v fll and v r . In this case muonic or tauonic asymmetries would be transformed 
through-oscillations into electronic asymmetry. This problem was analyzed recently 
in ref. 113 where it was shown that for large mixing angle solution for the solar 
neutrino anomaly and for (y^ — z^ T )-mixing deduced from the atmospheric neutrino 
anomaly the bounds on chemical potentials of all flavors are approximately 

16*1 < 0.1 (11) 

5 Neutrino Oscillations and BBN 

Effects of neutrino oscillations on BBN are very much different if only active neutri- 
nos are mixed or mixing is allowed between active and hypothetical sterile neutrinos. 
In the first case oscillations do not create any deviation from the standard BBN 
results if neutrinos are in thermal equilibrium with vanishing chemical potentials. 
Indeed, in this case oscillations would not lead to any modification of the standard 
distribution functions of neutrinos and abundances of light elements would remain 
the same as they were without oscillations. A noticeable effect would arise if neutri- 
nos are strongly degenerate and different lepton asymmetries would be redistributed 
by oscillations as discussed in the previous section. 

Physics is much more interesting if there are oscillations between active and 
sterile neutrinos. One evident effect is that oscillations would create additional 
neutrino species leading to N v > 3. Second, oscillations may distort spectrum of 
electronic neutrinos. The sign of the effect may be both positive or negative de- 
pending upon the form of spectral distortion. Third, in the case of MSW-resonance 
oscillations between v e and v s could be more efficient than oscillation between an- 
tineutrinos or vice versa. This would give rise to generation of lepton asymmetry 
in the sector of active neutrinos and, in particular, of electronic charge asymmetry 
which would have a strong impact on BBN. 

Excitation of additional, sterile, degrees of freedom at BBN by neutrino oscil- 
lations was considered in many papers starting from 1990 (a large list of references 
can be found in 111) . A recent bound on_the mixing parameters between v s and 24 
and v a and or v r respectively reads li-3: 

(6ml„JeV 2 ) sin 4 2^' = 3.16 • 10- 5 ( 5 »(T;; od )/10.75) 3 (AA,) 2 (12) 
(Sm^JeV 2 ) sin 4 W v <£ = 1.74 • 10" V(r;; od )/10.75) 3 (Ai\g 2 (13) 

where the number of relativistic degrees of freedom is taken at the temperature 
Tprod a ^ wmc h sterile neutrinos are effectively produced: 

T v p ; od = (12, 15) (S/y) 1 ^ (JmVeV 2 ) 1 ^ MeV (14) 

These bounds are valid only if AiV„ < 1 and spectral distortion of v e is neglected. 

The impact on BBN of the distortion of the spectrum of v e by oscillations was 
discussed e.g. in ref. Il3. According to this work an analytical fit to the bound on 
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the oscillation parameters that follows from the consideration of primordial 4 He 
can be written as 

Sm 2 (sin 2 26») 4 < 1.5 • 1CT 9 eV 2 , for Sm 2 < l(T 7 eV 2 . (15) 

The situation is much more complicated for a larger mass difference. In par- 
ticular, for the mass difference Sm 2 *~ (1 — 100) eV 2 and a small vacuum mixing 
angle, sin 2 29 < 10~ 3 , the resonance amplification of lepton asymmetry in the sec- 
tor of active neutrinos can take place £3. The effect was discussed in detail in the 
subsequent literature and now seems to be confirmed both numerically and analyt- 
ically. The impact of this phenomenon on BBN could be quite strong but no simple 
analytical results have been presented. An example of calculation of the effective 
neutrino number AN U induced by the generation of asymmetry in v e sector can be 
found in ref. Ell. 

6 Spatial Variation of Primordial Abundances 

An interesting phenomenon would arise if lepton asymmetries were inhomogeneous 
and large during BBN. In this case a large spatial variation of primordial abun- 
dances of light elements should take place. If the scale of variation is larger than 
the mixing scale (galactic) , then the observed abundances of light elements would 
be different in different regions of the universe. A mechanism of generation of large 
and inhomogeneous lepton asymmetry together with a small, and possibly homoge- 
neous baryon asymmetry was considered long ago in ref. E3 in the frameworks of 
Affleck and Dine lepto/baryogenesis scenario. Another model of creation of large 
and inhomogeneous lepton asymmetry by resonance neutrino oscillations (similar 
to the discussed in the previous section) was suggested recently in the paper ta. 

To avoid large density perturbations generated by varying chemical potentials 
it was assumed E3 that there exists symmetry with respect to permutation of elec- 
tron, muon and tauon asymmetries. In the simplest version of such model elec- 
tron asymmetry is small over 2 /3 of the sky but muon or tauon asymmetries are 
large (of order unity). Abundances of light elements are normal there, i.e. the 
mass fraction of primordial 4 He is Y p sa 0.25 and the deuterium-hydrogen ratio is 
D/H = 3 • 10~ 5 . Over 1/6 of the sky where electron asymmetry is large and nega- 
tive and other asymmetries are small the primordial abundances are high, Y p « 0.5 
and D/H = 10 • 10~ 5 , and over other 1/6 of the sky where electron asymmetry is 
large and positive the abundances are low, Y p « 0.12 and D/H = 1.5 • 10~ 5 . In 
more complicated versions of the model the probability distribution of the regions 
with normal, high, and low primordial abundances could be different. 

The characteristic scale of variation is not predicted by the model but from the 
upper limit on angular variation of the temperature of GMBR one may conclude 
that the scale should be larger than a few hundreds Mpc c3. Still some features in 
the angular spectrum of CMBR may be observable, in particular, diffusion damping 
slope at high I could be different in different directions in the sky. 

The model discussed above presents an example of cosmology with broken 
Copernicus Principle: the universe is energetically smooth but strongly chemically 
inhomogeneous. 
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7 Non-standard Neutrino Properties 



If neutrinos are massive (with Dirac mass) or possess right-current interactions then 
in addition to the usual left-handed species right-handed neutrino states might be 
present at BBN. Kinematical excitation of right-handed neutrino states by their 
Dirac mass was considered in many papers (see the review El). The latest and the 
most accurate work presents the limits: 



< 



130 keV, T QCD = 100 MeV 
120 keV, Tqcd = 200 MeV 



j 150 keV, T QCD = 100MeV . . 

- > 140 keV, Tqcd = 200 MeV. { ' 



if AiV„ < 0.3, while for AN U < 1.0, they are 



< 



310 keV, Tqcd = 100 MeV 
290 keV, Tqcd = 200 MeV 



f 370 keV, Tqcd — 100 MeV 
" T - \ 340 keV, Tqcd — 200 MeV. [ ' 

These limits are much stronger than laboratory limits for v T mass and comparable 
to the limit on mass. They are applicable if the neutrino life-time is longer than 
the characteristic time of nucleosynthesis (of course, because of Gerstein-Zeldovich 
limit E3 the life-time should be shorter than the universe age). On the other hand, 
an interpretation of neutrino anomalies in terms of neutrino oscillations demands 
very small mass difference and together with the laboratory bound, m„ e < 3 eV 
(see ref. E3), they lead to the masses of and v T much below the bounds ( fil||l7|) . 

If neutrinos, in addition to left-handed currents, are also coupled to right-handed 
ones, then the right-handed degrees of freedom could be present at BBN even if 
neutrinos are strictly massless. To avoid too many right-handed neutrinos the inter- 
actions with right-handed currents should be sufficiently weak and decouple at high 
temperatures such that the entropy dilution due to massive particle annihilation 
would diminish the number density of right-handed neutrinos down to a safe for 
BBN value. In other words, the mass of right-handed intermediate bosons should 
be sufficiently high. There is some disagreement in the literature about the the 
exact value of the lower bound on_their mass, but roughly speaking it should be 
larger than a few TeV (see review Id) . 

Another source of production of right-handed neutrinos could be their magnetic 
moment. There are two possible types of processes in the early universe where 
neutrino spin flip due to magnetic moment might take place: first, direct particle 
reactions, either quasi-elastic scattering, + vl — > + vr, and annihilation, 
e~ + e + — > vl,r + vr.l, or the plasmon decay, j p i — > vl,r + ^r,l- The second 
process is the classical spin rotation of neutrinos in large scale primordial magnetic 
fields that might have-existed in the early universe. The—former mechanism was 
first considered in ref. cJ, while the second one in refs. EBcB. 

Assuming that the extra number of neutrino species allowed by BBN is smaller 
than 1, the authors of the most recent paper E3 found the limit: 

H V < 2.9 • 10- 1( Vs (18) 
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for the case of reaction produced neutrinos. 

Potentially stronger limit can be obtained from consideration of neutrino spin- 
flip in cosmological magnetic fields at nucleosynthesis epoch. Unfortunately the 
result depends upon poorly known field strength in the early universe. With a 
reasonable assumption about the latter the upper bound on magnetic moment of 
neutrinos could be much stronger than ( |l8| ) , see discussion in the review El. 
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